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1 Introduction 

The i?-parity (Rp) symmetry plays an essential role in the construction of supersym- 
metric theories, such as the Minimal Supersymmetric extension of the Standard Model 
(MSSM) [1]. The conservation of Rp is closely related to the conservation of lepton (L) 
and baryon (B) numbers and the multiplicative quantum number associated to the Rp 
symmetry is defined by Rp = (_ i)3-B+l+25 j^j, ^ particle with spin S [2]. Standard model 
particles have even Rp, whereas the corresponding superpartners have odd Rp. The 
conservation of Rp guarantees that the spin-0 sfermions cannot be directly exchanged 
between standard fermions. It also implies that the sparticles {Rp = — 1) can only be 
produced in pairs, and that the decay of a sparticle leads to another sparticle, or an odd 
number of them. Therefore, it ensures the stability of the Lightest Supersymmetric Par- 
ticle (LSP). In the MSSM, the conservation of Rp is assumed: this is phenomenologically 
justified by proton decay constraints, and by the fact that a neutral LSP could be a good 
dark matter candidate. 

From a theoretical point of view, the conservation of Rp is not mandatory in super- 
symmetric extensions of the Standard Model (SM). Nevertheless, to be in agreement 
with the present experimental limit on proton lifetime, Rp violation can be introduced 
in MSSM either via the non-conservation of L or the non-conservation of B. One of the 
major consequences of the non-conservation of Rp is the allowed decay of the LSP into 
fermions; this modifies the signatures of supersymmetric particle production compared 
to the expected signatures in the case of Rp conservation. 

In this paper, searches for pair-produced supersymmetric particles in the hypothesis 
of Rp violation via one dominant sparticle-particle coupling are presented. The data 
recorded in 1999 and 2000 by the DELPHI experiment have been analyzed, and no signal 
of -Rp-violating decays was found in any of the channels. Previous results published by 
DELPHI on this subject can be found in references [3,4]. Similar searches performed by 
the other three LEP experiments have also shown no evidence for i?p- violating effects [5] . 

The paper is organized as follows. Section 2 is dedicated to the Rp violation phe- 
nomenology considered in the present search. The data samples and simulated sets are 
presented in Section 3. Section 4 is devoted to the description of the analyses, and in 
Section 5 the search results are given and interpreted in order to constrain the mass 
spectrum of SUSY particles. A brief summary is given in the last section. 

2 Rp non-conservation framework 

In the presence of Rp violation the superpotential [6] contains three trilinear terms, 
two violating L conservation, and one violating B conservation. We consider here the 
AjjfcLjLjEfc (non-conservation of L) and X'--fJ]iI)jI)k (non-conservation of B) terms ^, 
which couple the sleptons to the leptons and the squarks to the quarks, respectively. 
Since A^-fc = - Xjik and A^'^-^ = -A^'^^-, due to SU(2) and SU(3) symmetries, there are only 
9 Xijk and only 9 X'-^^ free couplings. In the present work, it is assumed that only one Xijk 
or X'-jj, is dominant at a time. In the following, searches assuming i?p-violation via one 
dominant AyfcLjLjEfc term are referred as "LLE", and those via one X'-jiJJiT)jt)k term as 
"UDD". Searches assuming i?p-violation via one A^ ^.LjQjDfc term (non-conservation of 
L) were not performed in DELPHI for data collected in 1999 and 2000. 

^i,j, k arc generation indices, L denotes the lepton doublet superfields, E {U, D) denote the lepton (up and down quark) 
singlet superfields, \ijk and A" ,. are Yukawa couplings. 



In the pair-production of supersymmetric particles studied here, Rp is not conserved 
in the decay of the sparticles, but is conserved at the production vertex. The production 
cross-sections behave as in the MSSM with Rp conservation (see section 2.3). 

2.1 i?p-violating decays of sparticles via LLE or tJDD terms 

Two types of supersymmetric particle decays are considered: direct decay and indirect 
decay. 

2.1.1 Direct decays 

Rp violation allows the direct decay of a sfermion into two conventional fermions 
(Fig. 1-a, b), or the direct decay of a neutralino or a chargino into a fermion and a 
virtual sfermion which then decays into two conventional fermions (Fig. 1-c). A direct 
decay is the only possibility for the LSP. 

Decays through LLE terms 

Sleptons are coupled to leptons through the AjjfcLjLjEfc term. In four-component Dirac 
notation, the LLE Yukawa interaction terms are^ [7]: 

\jk (vidkdjL + ijdkRViL + ^kRO^iLY^jL -i^ jj + h.c. 

Considering the above expression, it can be deduced that the i?p-violating decay of a 
sfermion is possible only with specific indices i,j,k of the coupling which is considered to 
be dominant. The possible sparticle decays with such a dominant Xijk coupling are listed 
below. 

. The sneutrino direct decay gives two charged leptons: via Xijk only the Ui and Uj 

are allowed to decay directly: z/j -^ ^fL^^R ^^(^ ^j ~^ ^fh^'kR respectively. 
. The charged slepton direct decay gives one neutrino and one charged lepton (the 
lepton flavour may be different from the slepton one). Among the supersymmetric 
partners of the right-handed leptons, only the one belonging to the k^^ generation 
can decay directly: £^^ -^ ^iL^^l > ^Tl^jl- For the supersymmetric partners of the 

left-handed leptons, the allowed direct decays are: £j^ — > i^jl^^r ^^^ ^jL ~^ ^iL^m- 
. The neutralino decays via a virtual slepton and a lepton, and subsequently gives 
three-lepton final states (two charged leptons and one neutrino): 

. The chargino decays via a virtual slepton and gives either three charged leptons, or 
two neutrinos and one charged lepton: 

Decays through UDD terms 

The squarks are coupled to the quarks through the \'-jiJJiT)jT)k term. The decays allowed 
via this term can be inferred by considering the Lagrangian for the trilinear Yukawa 
interactions written in expanded notation: 

a:;., [w jd^Tdl +Wd*Wr + n* w Wr) + h.c. 

From this Lagrangian, we can derive the following rules: 

^Hcrc u (v) refer to neutrino (sneutrino) fields, t. (I) refer to charged lepton (slepton) fields, i,j,k are generation indiees 
and the superfix c refers to a charge-conjugate field. 



. The direct decays of squarks into two quarks are given by: Ui^R -^ (ij,R(ik,R, 

dj,R — > ^i,Rdk,R and dk,R — > ^i,Rdj^R. 
. The neutrahno decays via a virtual squark and a quark and subsequently gives a 

three-quark final state: 

X°^ Ujdjdk , Ujdjdk 
. Chargino decay is similar to the neutralino one, and then gives also a three quarks 

final state: x^— > UjdjUk , UjUjdfc , didjdk 

2.1.2 Indirect decays 

Indirect decays are cascade decays through i?p-conserving vertices to on-shell super- 
symmetric particles, down to the lightest supersymmetric particle, which then decays 
via one LLE or UDD term (Fig. 2). A typical example is the i?p-conserving decay 
Xi^' Xi+ W*+ (see Fig. 2-e) and the subsequent decay of x? through the i?p-violating 
couplings. The indirect decay mode usually dominates when there is enough phase space 
available in the decay between "mother" and "daughter" sparticles. For example, when 
the difference of masses between these two sparticles is larger than 5-10 GeV/c^. Re- 
gions of the parameter space where there is a "dynamic" suppression of the i?p-conserving 
modes also exist. In this case, even if the sparticle is not the LSP, it decays through an 
i?p-violating mode. For example, if the field component of the two lightest neutralinos is 
mainly the photino, then the decay X2 ~^ Xi Z* is suppressed. 

The sfermion indirect^ decay studied here is the decay through the lightest neutralino 
considered as the LSP (f^ f Xi); followed by the -Rp- violating decay of the LSP. With 
LLE, the indirect decay of a sneutrino (charged slepton) through a neutralino and a 
neutrino (charged lepton) leads to two charged leptons and two neutrinos (three charged 
leptons and one neutrino). The squark decay into a quark and a gaugino leads to one 
quark and three leptons. With IJDD, the indirect decay of a squark (slepton) leads to 
four quarks (three quarks and one lepton). 

2.2 i?p-violating coupling upper limits and LSP lifetime 

Upper limits on the Xijk and A^'-^ couplings can be derived mainly from indirect searches 
of i?p- violating effects [7,8], assuming that only one coupling is dominant at a time. They 
are dependent on the sfermion mass, and usually given for mj = 100 GeV/c^. The upper 
bounds on Xij^ are obtained from charged-current universality, lepton universality, t'^— e 
scattering, forward-backward asymmetry in e"'"e~ collisions, and bounds on z/g-Majorana 
mass. Most present indirect limits are in the range of 10~^ to 10"^; the most stringent 
upper limit is given for A133 (~ 6 ■ 10^^). Upper limits on A^'-^ couplings come from 
experimental measurements of double nucleon decays for A'/^2 (lO^^); n — n oscillations 
for A'/i3 (10-5) and of Rt = Thad{^°)/Ti{Z^) in e+e^ collisions for A'3'12, A'g'ig, A^j'as (0-43). 
The upper limits on the other A^' ;, couplings are obtained from the requirement of per- 
turbative unification at the Grand Unified Theory (GUT) scale of 10^^ GeV. This gives 
a limit of 1.25. 

In the present searches, the LSP lifetime was a crucial parameter since the analyses 
were valid only if the i?p-violating decays were close to the production vertex, which 
means a LSP fiight path shorter than a few centimetres. 



The LSP mean decay length is given by [9,10]: 



100 GeV/cV \^ m^ J A'^ 

if the neutrahno or the chargino is the LSP with f3'j = P^/m^^ and with A = X^k or 
A = \^\'-ji^. Considering the upper hmits on the couphngs described above and according 
to equation (1), the analyses are not sensitive to a light neutralino (m^ < 15 GeV/c^), 
due to the terms m^^"^ and (/37). Moreover, when studying neutralino decays, for the 
typical masses considered in the present study, the analyses are sensitivite to i?p-violating 
couplings greater than 10""^ to 10~^, where the i?p-violating decay has a negligible decay 
length. For much lower values of the coupling strength, the LSP escapes the tracking 
devices before decaying and the results of the searches performed under the assumption 
of Rp conservation are recovered [11]. Between these two extreme cases, the LSP decay 
produces a displaced vertex topology ^. 

2.3 Pair-production of supersymmetric particles 

Pair-production of supersymmetric particles in the MSSM assuming Rp violation is 
identical to pair-production in the case of Rp conservation, since the trilinear couplings 
are not present at the production vertex. The production of single supersymmetric parti- 
cles via trilinear couplings has been studied in other searches performed by the DELPHI 
collaboration [13]. 

In the constrained MSSM scheme [1] considered in the present searches, the mass 
spectrum of neutralinos and charginos is determined by three parameters, with the as- 
sumption that both the gaugino and the sfermion masses are unified at the GUT scale. 
The relevant parameters are then: M2, the SU(2) gaugino mass at the electroweak scale 
(it is assumed that Mi = |tan^^vvM2), mg, the common sfermion mass at the GUT scale, 
/i, the mass-mixing term of the Higgs doublets at the electroweak scale and tan/5, the 
ratio of the vacuum expectation values of the two Higgs doublets. It is assumed that the 
running of the A^fc and A^'^^ couplings from the GUT to the electroweak scales does not 
have a significant effect on the "running" of the gaugino and sfermion masses. 

The charginos are produced in pairs in the s-channel via 7 or Z exchange, or in the 
t-channel via z/g exchange if the charginos have a gaugino component; the neutralinos are 
produced in pairs via s-channel Z exchange provided they have a higgsino component, or 
via t-channel e exchange if they have a gaugino component. The t-channel contribution 
is suppressed when the slepton masses (depending on mg) are high enough. When the e 
mass is sufficiently small (less than 100 GeV/c^), neutralino production can be enhanced, 
because of the t-channel contribution. On the contrary, if the z/g mass is in the same 
range, the chargino cross-section can decrease due to destructive interference between 
the s- and t-channel amplitudes. 

The pair-production cross-section of sfermions mainly depends on the sfermion masses. 
The e and u^ cross-sections are also very sensitive to the neutralino and chargino compo- 
sitions (which are function of fi, M2 and tan/5) via the t-channel exchange. The sfermion 
mass-eigenstates, fi and f2 (where f is a quark or lepton and fi is lighter than f2), are 

■^This particuliar topology, not considered in the present searches, has been studied in other searches performed by the 
DELPHI collaboration [12]. 



obtained from the two supersymmetric scalar partners fj, and fji of the corresponding left 
and right-handed fermion [14,15]: 

fi = h cos$mix +Jr sin<l>mix 
h = -h sin$mix + fi? cos<l> 



mix 



where $mix is the mixing angle with < $mix ^ ^r. The supersymmetric partner of the 
left-handed fermions are likely to be heavier than their right-handed counterparts. The 
fL^fn mixing is related to the off-diagonal terms of the scalar squared-mass matrix. It 
is proportional to the fermion mass, and is small compared to the diagonal terms, with 
the possible exception of the third family sfermion [16]. The lighter stop, ti, is then 
probably the lightest squark. This is not only due to the mixing effect but also to the 
effect of the large Yukawa coupling of the top; both tend to decrease the mass of ii [17]. 
Similarly the lightest charged slepton is probably the fi. For small values of tan/?, fi is 
predominantly a tr, and it is not so much lighter than Ci and jli. In the present slepton 
search, a no-mixing scenario is assumed. In the third squark generation searches two 
left-right mixing angle cases have been considered. The first one with mixing angle equal 
to zero and the second one with the mixing angle <l>mix = 56° ($mix = 68°) corresponding 
to the minimum production cross-section of the stop (sbottom) via Z exchange [18]. 

3 Data and generated samples 

3.1 Data samples 

The data recorded in 1999 and 2000 by the DELPHI experiment at centre-of-mass 
energies from y/s = 192 GeV to 208.8 GeV, correspond to a total integrated luminosity 
of around 450 pb^^. The DELPHI detector has been described elsewhere [19]. An 
integrated luminosity of 386 pb^^ (Table 1) has been analysed, corresponding to high 
quality data, with the tracking detectors and the electromagnetic calorimeters in good 
working conditions. At the end of the data taking period in 2000, one sector (among 
twelve) of the Time Projection Chamber (TPC) failed beyond repair. This required 
modifications in the data treatment (pattern recognition), and a specific simulation of 
the detector with one TPC sector off has been performed. An integrated luminosity of 
51.8 pb^^ recorded with one TPC sector off have been analysed. 

3.2 Event generators 

To evaluate background contaminations, different contributions coming from the SM 
processes were considered. The SM events were produced by the following generators: 
. 77 events: BDK [20] for 77 -^ i^£~ processes, including radiative corrections for the 

e"'"e~/i"'"/i~ and e"'"e~r+r~ final states, and TWOGAM for 77 — > hadron processes. 
. two-fermion processes: BHWIDE [21] for Bhabha scattering (e"'"e~ -^ e"'"e~(7)), 
KORALZ [22] for e+e" -^ AiV~(7) and for e+e" -* r+r~(7) and PYTHIA 6.143 [23] 
for e+e~ — ^ qq(7) events. 
. four-fermion processes: EXCALIBUR [24] and GRC4F [25] for all types of four-fermion 
processes: non resonant (fff'f'), singly resonant (Zff, Wff') and doubly resonant (ZZ, 
WW) (PYTHIA was used also for cross-checks on the final results). 
Signal events were generated for all analyses with the SUSYGEN 3 . 00 program [26] . 

All generated background and signal events were passed through the full DELPHI 
simulation and reconstruction chain [19] and then processed in the same way as the real 



data. To treat the data taken with one sector of the TPC off, special background and 
signal event samples were generated, and the same treatment applied to them as to the 
real data. 

3.3 Signal samples 

Choice of the i?p-violating couplings 

Among the nine Xijk couplings, A 122 (which leads to several muons in the final states) and 
A133 (which leads to several taus in the final states) have been chosen for most of the sig- 
nal generation. Their values were set for m^ = 100 GeV/c^ at 0.04 and 0.003 respectively, 
below their upper bound derived from indirect searches of i?p-violating effects. Any value 
between 10~^ and 10~^ would not change the neutralino decay topologies. Simulations 
with other couplings have been also performed in order to check that the analyses devel- 
oped for A122 or A133 were able to select the corresponding signal with an equal or better 
efficiency. 

For the generation of all IJDD signals, a A212 coupling of strength 0.1 (for mq = 
100 GeV/c^) was used. Any value between 10~^ and 0.5 would not change the neutralino 
decay topologies. 

Searches for decays through specific A^'^-^ couplings, leading to the production of one 
or several b quarks, can use b-tagging techniques to reach higher sensitivities, but at the 
cost of losing generality. 

Generated signal sets 

Two different procedures were applied to the signal generation for gaugino pair-production 
and subsequent decays through either LLE or IJDD terms in order to cover the MSSM 
parameter space. 

For the LLE term, the x° and Xk pair-production processes were considered for differ- 
ent values of tan/9 (from 1 to 30), mo (between 90 GeV/c^ and 500 GeV/c^), fi (between 
-200 GeV/c^ and 200 GeV/c^) and M2 (between 5 and 400 GeV/c^), for centre-of-mass 
energies of 200 and 206 GeV. For the UDD term, pair-production of neutralinos was 
generated for several masses. The simulated masses started from 10 GeV/c^ and were in- 
creased in steps of 10 GeV/c^, as long as the mass of the chargino remained kinematically 
accessible. Masses corresponding to the kinematic limit were also simulated. To generate 
chargino pairs, the mass of the chargino was varied from 45 GeV/c^ to 95 GeV/c^ with 
a 10 GeV/c^ step. Chargino masses were also simulated at the kinematic limit. The 
neutralino mass was varied from 10 GeV/c^ to a mass difference with the chargino of 
5 GeV/c^ with a 10 GeV/c^ step. For each mass pair, a set of the variables /x, M2 and 
tan/5 was found for the chosen simulation. 

Sfermion indirect decay signals were simulated at different masses with steps of 
10 GeV/c^ at centre-of-mass energies of 200 and 206 GeV, with tan/3 and n fixed at 
1.5 and -200 GeV/c^ respectively. M2 was used to fix the neutralino mass at the required 
value. The points were simulated from 45 to 100 GeV/c^ for the sfermion masses and 
from 15 to 95 GeV/c^ for the Xi masses up to a mass difference between the sfermion 
and the LSP of 5 GeV/c^. 

Among the sfermions, only the sneutrino direct decay via LLE terms was studied. 
Specific signal sets have been produced with Br(z/ -^ I'^Cr) = 100%. The processes 



z/gZ/e^ 4/i (A122), i^ei^c-^ 4:T (A133), u^u^^ At (A233) and v^v^^ 2e2r (A133) have been 
generated for different values of the sneutrino mass up to 98 GeV/c^, with tan/3 and /x 
fixed at 1.5 and -200 GeV/c^ respectively. In order to check that all final states from 
vv decay were covered, signals obtained for other Xijk couplings and for sneutrino masses 
around 90 GeV/c^ were also generated. 

4 Description of the analyses 

The analyses covering the decay of pair-produced sparticles were designed to cover 
multi-lepton final states for LLE coupling and multi-jet final states for UDD coupling. 
Different preselections were applied, one for the multi-lepton channels and one for the 
multi-jet channels. In each case, dedicated analyses were necessary to take into account 
the specific characteristics of the sparticle decay. The multi-jet analyses required a specific 
treatment based on neural network techniques. 

The sensitivity of the searches for sparticle indirect decays depended on the mass 
difference (AM) between the sparticle being searched for and the LSP. The analyses were 
designed to be efficient for AM > 5 GeV/c^. The multi-jet analyses required different 
signal selection optimisations to cover efficiently all AM regions; therefore they were 
divided into windows according to the value of AM. 

No excess in the data appeared in these searches, therefore a working point optimiza- 
tion on the selection criteria was performed minimizing the expected excluded cross- 
section as a function of the average signal efficiency. 

4.1 Description of the final states 
4.1.1 Decays via LLE 

Direct and indirect decays of gauginos, direct and indirect decays of sneutrinos and 
indirect decays of charged sleptons and squarks were studied. 

The direct decay of a pair of lightest neutralinos leads to two neutrinos and four 
charged leptons. For an indirect decay of chargino or heavier neutralino pairs the final 
state may contain some jets and/or leptons in addition to the four leptons and the missing 
energy from the decay of the LSP. The direct decay of a sneutrino pair gives final states 
with four charged leptons, in which the leptons can be of two different flavours. The 
direct decay of a charged slepton pair gives flnal states with two charged leptons, in 
which the leptons can be of two different flavours, and missing energy. This flnal state 
has not been covered by the present analyses. In the indirect decay of any sfermion pair, 
the flnal states are composed of two fermions plus the decay products of the neutralinos. 

Compared to other couplings the highest efficiencies and background reduction were 
obtained in analyses performed on the signal with a dominant A122 coupling. For analyses 
dedicated to a A133 coupling, due to the presence of several taus in the decay channels, the 
efficiencies and the rejection power were low. For final states produced by other A^-fc, the 
detection efficiencies lay between these two limiting cases. Therefore conservative limits 
can be derived by considering the results of the analyses performed assuming a dominant 
A133 coupling, and only these analyses will be described in section 4.3. 

The decay of pair-produced sparticles via a A133 coupling leads to different types of 
final states, depending on the produced sparticles. The XiXi decay via A133 leads to 
2r + £ + £'-|-^ , where £, £' = e or r, and Ip means missing energy. In addition, jets and/or 
leptons from the W or Z decays show up in the final state from Xi'x? and xtXi indirect 



decays. The indirect decay of a slepton pair gives 2r + £ + (^'+Ifi , with two additional 
charged leptons (same flavour, opposite charge) in the case of charged sleptons, and 
additional missing energy in the case of sneutrinos. A 4r final state is produced by the 
direct decay of z/q^c via A133. The direct decay of z/^^r gives 2e2r and then, there is finally 
less missing energy comming from the taus decay than in the previous cases. The indirect 
decay of squarks adds exactly two jets to the 2r + £ + (i'+Ip . 

Four analyses have been performed to search for all these topologies. They are sum- 
marized in the first part of Table 2. 

4.1.2 Decays via UDD 

Direct and indirect decays of gauginos, and indirect decays of charginos, charged slep- 
tons and squarks were studied. 

For each indirect decay of a chargino, squark or slepton pair there are at least six 
quarks in the final state. Therefore the most important feature of these decays is the 
number of quarks produced, which can be up to ten for the indirect decay of two charginos 
with the hadronic decays of the W bosons. The indirect decay channel presents the only 
possibility for the sleptons to decay through a IJDD term. In this case, two leptons are 
produced in the Rp conserving decay of the slepton pair, and they add to the six jets 
coming from the decay of the two neutralinos: a 6 jets + 2£, £ = e, /i final state is the 
signature of these signals. The indirect decay of a stop or sbottom pair produces 8 jets 
in the final state. Two b quarks are produced from the sbottom decay. The analysis of 
the different decay channels was organized on the basis of the number of hadronic jets in 
the final state (see Table 2). 

4.2 Analysis tools and techniques 

4.2.1 Lepton identification 

The identification of a muon or an electron, used in all Xijk analyses and several Af •;. 
ones, was based on standard DELPHI algorithms [19]. The identification could be "tight" 
if it was an unambiguous one, or "loose" otherwise. In the multi-lepton analyses described 
in this section a particle was considered as a well identified electron if it satisfied the tight 
conditions from the DELPHI electron identification algorithm, its momentum was greater 
than 8 GeV/c and there was no other charged particle in a cone of half-angle 2° around 
it. A particle was considered as a well identified muon if its momentum was greater than 
5 GeV/c and it was tagged as a tight muon candidate by the DELPHI algorithm. 



4.2.2 Jet reconstruction algorithms 

Two different jet reconstruction algorithms have been used. The DURHAM algorithm [27] 
was used for the multi-lepton (LLE coupling) analyses, where jets were expected from r 
or W boson decays. In case of multi-jet analyses, the CAMBRIDGE clustering algorithm [28] 
implemented in the CKERN package [29] was used. 

The CAMBRIDGE algorithm was introduced to select soft jets, coming from quark-jets 
with gluon emission. The specific procedure of clusterization which extracts soft jets 
from the list of objects to be clustered, was particularly interesting for multi-jet analyses, 
where the jets (more than six) may not be well separated in momentum space. For 
each event, the two algorithms provided all possible configurations of jets between 2 and 



10. They have the same definition of jcut distance, but mainly differ in the iterative 
procedure of clustering. In this paper, the transition value of the ycut in the DURHAM or 
CAMBRIDGE algorithm at which the event changes from a clustering with n jets, called 
n-jet configuration, to a clustering with (n— l)-jets, is denoted jnn-i- In other words, the 
jnn-i value is the ycut value for which the number of particle clusters fiips from n to n—1 
for increasing y^ut distances. For example, the y43 value of one event is the highest value 
of ycut to obtain 4 separated clusters of particles. 

4.2.3 Neural networks 

A neural network method was applied in order to distinguish signals from SM back- 
ground events for all multi-jet analyses. The trainings of the neural networks were done 
in the standard back-propagation manner with one hidden layer on samples of simulated 
background (qq and four-fermion) and signal events. A feed-forward algorithm has been 
implemented to compute from the input discriminating variables a single discriminant 
variable (signal output node) which was used first to validate the training with different 
signal and background samples and then to select the final number of candidate events for 
each analysis. The exact configuration and the input discriminating variables of each neu- 
ral network depended on the search channel. The working point on the signal ouput node 
value has been chosen to minimize the expected excluded cross-section at 95% confidence 
level (CL) when there is no signal. 

4.3 Multi-lepton LLE channels 

4.3.1 Preselection 

The selections were based on the criteria already presented in [3], using mainly miss- 
ing momentum, lepton identification and kinematic properties. The preselection require- 
ments were: 

. more than three charged particles and at least one of them with a polar angle between 
40° and 140°; 

. at least one identified lepton (e or fj,); 

. a total energy greater than O.l-^s; 

. a missing momentum component transverse to the beam (pt) greater than 5 GeV/c; 

. a polar angle of the missing momentum (^miss) between 27° and 153°; 

. a thrust axis not close to the beam pipe, viz. | cos^tij less than 0.9; 

. an acollinearity*^ greater than 2°, and greater than 7° for events with a charged 
particle multiplicity greater than 6. 

The preselection was efficient in suppressing 99.9% of the backgrounds coming from 
Bhabha scattering and two-photon processes while removing 97% of the ff(7) contribu- 
tion. The preselection also reduced the four-fermion contamination by 75%. After this 
preselection stage, 2310 events (1220 for data at centre-of-mass energies between 192 and 
202 GeV, and 1090 for those collected above 202 GeV) were selected to be compared to 
2254 ± 6 expected from the background sources (1189 ± 5 at centre-of-mass energies 
between 192 and 202 GeV, 1065 ± 4 above 202 GeV). The corresponding efficiencies for 
the large majority of LLE signals lay between 60% and 80%. The distributions of several 
event variables at the hadronic preselection stage are shown in Figure 3. 

The above requirements had to be slightly modified for the stop analysis (see Section 
4.3.5), in particular to take into account the fact that the final state always contains two 

"^The acoUinearity is computed between the two vectors corresponding to the sum of the particle momenta in each 
hemisphere of the event. The two hemispheres are defined by the plane orthogonal to the thrust axis. 
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jets. A minimum multiplicity of eight charged particles was required. No requirement 
was applied on the thrust axis. On the other hand, a stronger cut was applied on the 
polar angle of the missing momentum (30° < ^miss < 150°). After this preselection stage, 
2197 events were selected to be compared to 2208 ± 6 expected from the background 
sources. 

4.3.2 Gaugino search 

The gaugino analysis was designed to cover the 2r + n£ + mj+ ^ (n > 2, m > 0) final 
states, from direct or indirect decays of gauginos, and to be efficient for both low and 
high multiplicity cases. 

The thrust had to be less than 0.9 and a lower limit on the missing energy was applied 
viz. Emiss greater than O.S-a/s. The number of neutral (charged) particles had to be less 
than 20 (25) and the polar angle of at least one lepton had to be between 40° and 140°. 

The events were then separated in to two classes, according to their charged particle 
multiplicity. 

• For events with a charged particle multiplicity from four to six, (mainly for neutralino 
direct decay topologies), the following criteria were applied: 

. the energy in a cone of 30° around the beam axis was required to be less than 

50% of the total visible energy; 
. the energy of the most energetic lepton (e or n) had to be between 2 and 70 GeV; 
. there should be no other charged particle in a cone of half angle of 20° (6°) 

around any identified lepton for a charged particle multiplicity equal to 4 (5 or 

6). 

• For events with a charged particle multiplicity greater than six, the previous criteria 
became: 

. the energy in a cone of 30° around the beam axis was required to be less than 

40% of the total visible energy; 
. the energy of the most energetic lepton had to be between 5 and 60 GeV, 
. if there was only one identified lepton; there should be no other charged particle 

in a cone with a half angle of 6° around it; if there was more than one identified 

lepton there should be no other charged particle in a cone with a half angle of 

10° around at least two of them; 
. at least one electron (loose identification) was required. 

These criteria removed 95% of ff7, ZZ and W"'"W~ events. 

A selection based on the jet characteristics and topologies was then applied, depending 
on the charged particle multiplicity, as mentioned above. First, constraints were imposed 
on the y32 and J43 values: they had to be greater than 0.002 and 0.0001 respectively for 
events with low charged particle multiplicity, and greater than 0.016 and 0.005 respec- 
tively for events whose charged particle multiplicity was above 7; these criteria eliminated 
99% of the remaining ff7 contribution. In events with more than six charged particles, 
it was required that at least one jet had no more than two charged particles. In four or 
five-jet configurations, a minimum number of 4 jets with at least one charged particle was 
required. For a four-jet topology, a cut was applied on the value of E^j^^ ■ 6'jnin where E:^;^ 
is the energy of the least energetic jet, and ^min is the minimum di-jet angle (Fig. 4-a). It 
had to be greater than 1 GeV-rad for events with low charged particle multiplicity, and 
greater then 5 GeV-rad for events with a charged particle multiplicity above 7. These 
requirements reduced the background from 4-fermion processes. 
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4.3.3 Slepton search 

A slepton analysis, aimed to search for 2r + 2£+ Ifi (+ 2€) was performed in order to 
studythe following three channels: 

. z/cZ/c -^ 4r. 
These final states have several taus, and most of them have missing energy. After the 
tau decays, all channels present a large amount of missing energy. The criteria used 
to eliminate almost all remaining ff7 events and most of the 4-fermion events were the 
following: 

. the missing energy had to be greater than 0.3-y^; 

. the energy in a cone of half-angle of 30° around the beam axis was required to be 
less than 40% of the total visible energy; 

. the number of charged (neutral) particles had to be less than 8 (10); 

. the energy of the most energetic lepton had to be between 2 and 70 GeV; 

. at least one lepton should have a polar angle between 40° and 140°; 

. there should be no other charged particle in a cone with a half-angle of 6° around 
at least one lepton; 

. y32 and J43, computed with the DURHAM algorithm, (Fig. 4-b) had to be greater than 
2 ■ 10"^ and 4 • 10~^ respectively; 

. in a four-jet topology a minimum angle of 20° between any pair of jets was required. 

4.3.4 Ur search 

An analysis searching for 2e2r final states produced in the direct decay of z/^z/^ was 
performed. Compared to the selection described in 4.3.3, the most important change 
was the suppression of the criterion on the missing energy, and the introduction of the 
requirement of having at least one well identified electron: 

. the energy in a cone of 30° around the beam axis was required to be less than 50% 
of the total visible energy; 

. the number of charged (neutral) particles had to be less than 7 (10); 

. there should be at least one electron; 

. the energy of the most energetic lepton had to be between 25 and 80 GeV; 

. at least one lepton should have a polar angle between 40° and 140°; 

. there should be no other charged particle in a cone of half-angle of 6° around at least 
one lepton; 

• y32 (Fig. 4-c), computed with the DURHAM algorithm, had to be greater than 2 ■ 10~^; 

. in a four-jet topology, a minimum angle of 20° between any pair of jets was required. 
This selection removed ff7 background and most of the 4-fermion events. 

4.3.5 Stop search 

In stop pair-production, each of the stops decays into a charm quark and a neutralino. 
The subsequent i?p-violating decay via the A133 coupling of the neutralino into leptons, 
the final state: 2r + 2£ + ^ +2j,i=e or r. 

After the preselection described in Section 4.3.1, the criteria used to select the final 
states of the stop pair indirect decay were: 

. the missing energy greater than 0.3■^/s^, 

. the charged and neutral particle multiplicities below 25 and 20 respectively; 
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. the polar angle of at least one lepton between 40° and 140°; 

. the energy of the most energetic identified lepton between 5 and 50 GeV; 

. no other charged particle in a cone of half-angle of 6° around at least one lepton; 

. at least one well identified electron (Fig. 4-d); 

. the y32 and 754 values constrained to be less than 0.016 and 10"'^ respectively; 

. in the four-jet configuration, at least one jet with less than three charged particles. 

4.4 Multi-jet tJDD channels 

For all UDD channels the main SM backgrounds come from four-fermion processes 
except for the low neutralino mass channel where the hadronic Z decay is the dominant 
background. The following IJDD analyses were based on neural network techniques since 
the optimisation of the signal selection over the four-fermion background was performed 
on topological variables, such as jet resolution parameters, which are extensively corre- 
lated. 

4.4.1 Preselection 

The multi-jet IJDD signals have final states with a large hadronic activity, independent 
of the produced sparticles. Therefore a general hadronic preselection was performed with 
the aim of a high efficiency for the signal (especially for the gaugino analysis) and at the 
same time a good rejection of low particle multiplicity hadronic background events: 

. the number of charged particles had to be greater than 15; 

. the total energy was required to be greater than O.G-a/s; 

. the energy associated to charged particles was required to be greater than 0.3-y/s; 

. the effective centre-of-mass energy ^ had to be greater than 150 GeV; 

. the discriminating variable da = amin ■ Emjn — 0.5 ■ /3min ■ Emax/Emin (where the 0.5 
energy factor is in GeV) had to be greater than -10 GeV.rad ^; 

. the minimum jet invariant mass had to be greater than 500 MeV/c^ when forcing 
the event into four jets; 

. the ln(y32) had to be greater than -6.9; 

. the ln(y43) had to be greater than -8. 

After the hadronic preselection, the main remaining background events were the four- 
fermion events and the qq7 events with hard gluon radiation. We observed 3844 events in 
the data with 3869 ± 4 events expected from background processes for the year 2000 (4180 
events in the data to be compared to 4096 ± 7 events in the simulation for the year 1999). 
Examples of the distributions of several event variables at the hadronic preselection level 
are shown in Figure 5. 

The efficiencies for the UDD signals varied from 60% to 99% depending on the simu- 
lated masses and AM. This preselection (sometimes with slight modifications, described 
in the following) was used for all the UDD analyses. 

These requirements had to be slightly modified to be better optimised for slepton 
searches. The discriminating variable do, was not used in the preselection, and the effective 
centre-of-mass energy was required to be above 0.6■^/s, because a tighter cut was set on 
y43 and a cut on y54 was applied. 

Additional criteria to the basic hadronic preselection have been applied before the 
optimal neural network selection: 

■'the effective centre-of-mass energy is the centrc-of-mass energy after the emission of one or more photons from the 
initial state. 

^c*min is the minimum angle between two jets, /^min is the minimum angle between the most energetic jet and any other, 
Emin (Eniax) is the minimum (maximum) jet energy from the four-jet topology of the event. 
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. the charged particle multiphcity had to be greater than 16; 

. the total energy was required to be greater than 0.6-y/s; 

. the energy of charged particles was required to be greater than 0.3-v^; 

. the ln(y43) had to be greater than -7; 

. the ln(y54) had to be greater than -8; 

. the thrust had to be lower than 0.94; 

. the maximum di-jet mass in a four-jet configuration had to be greater than 
10 GeV/c2. 
This selection was applied for the UDD sleptons ans squark analyses. After this selection 
the remaining number of events for all energies was 4245 for the data and 4378 ± 8 for 
the expected background from SM processes. 

The distributions of the variables ln(y54) and ln(y65) at this stage of the preselection level 
are shown in Figure 6. 



4.4.2 Neutralino search 

The analysis described here was mainly designed to search for neutralino direct decays; 
it was also efficient in the search for chargino direct decays. The six-jet analyses were 
based on three different neural networks for the optimization of the background and sig- 
nal discrimination. The neural network method used has been presented in Section 4.2.3. 

Events with low gaugino mass have a large boost and look like di-jet events. On the 
contrary, events with heavy gauginos are almost spherical with six well separated jets. 
Therefore, we distinguished 3 mass windows to increase the sensitivity of each signal 
configuration: 

. low mass window Nl: 10 < m^ < 45 GeV/c^; 

. medium mass window N2: 45 < m^ < 75 GeV/c^; 

. high mass window N3: m^ > 75 GeV/c^. 

A mass reconstruction was performed using a method depending on the mass window. 
For the Nl analysis, the events were forced into two jets and the average of the two- 
jet masses was computed. For the other analyses, the events were forced into six jets 
and criteria on di-jet angles were applied to choose the optimum three-jet combinations 
corresponding to the decays of two neutralinos with the same mass. The minimum 
and maximum angles between the jets belonging to the same three-jet cluster (same 
neutralino) had to be in the intervals [20°, 80°] and [50°, 165°] for the medium mass 
window N2 ([40°, 110°] and [100°, 175°] for the high mass window N3). If more than one 
combination was selected, the combination with the minimum difference between the two 
energies of the three-jet clusters was chosen to compute the neutralino mass. 

Three neural networks were used, one for each mass window, with the following vari- 
ables as inputs: 

. the thrust; 

. distww= \/ ^^-2 + ^^ ^:7 — ^^, where Mi and M2 are the di-jet masses of the 
V ""- '^+ 

jet combination which minimized this variable (after forcing the event into 4 jets); we 

took Mw = 80.4 GeV/c2 for the W mass, ct_ = 9.5 GeV/c^ and ct+ = 4.8 GeV/c^ 

for the mass resolutions of the difference and the sum of the reconstructed di-jet 

masses respectively; this variable is peaked at for WW events, allowing a good 

discrimination against this background; 
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. the energy of the least energetic jet muhiphed by the minimum di-jet angle in four 

and five-jet configurations; 
. the difference between the energies of the two combinations of three jets, after the 

mass reconstruction; 
. the reconstructed neutralino mass; 

• Yn n-i with n=3 to 10. 

It was observed that the modelling of the gluon emission was unable to describe the 
event distributions of y„ „„i correctly for n greater than 5. To take into account this 
imperfect description, corrections were applied to the background distributions of the 
Yn n-1 variables (for n between 6 and 10) [30]. 

4.4.3 Chargino search 

To take into account the effect of the mass difference between chargino and neutralino, 
AM, on the topology of the event, the ten-jet analysis was divided into two windows: 

. low AM window CI: 5< AM< 10 GeV/c^; 

. high AM window C2: AM> 10 GeV/c^ 

Two neural networks were trained with the following discriminating variables as in- 
puts: 

. the thrust; 

. the variable distww described above; 

. the energy of the least energetic jet multiplied by the minimum di-jet angle in four 
and five-jet configurations; 

• Jn n-i with n=3 to 10 (for n=6 to 10, the corrected y„ „_i were used); for the CI 
analysis, the variables js7, Jqs and yiog were not used. 

4.4.4 Slepton search 

Three mass windows were defined to take into account the mass difference, AM, be- 
tween the sfermion and the neutralino considered as the LSP: 

. low AM window 1: 5 < AM < 10 GeV/c^ with m^o > 55 GeV/c^; 

. high AM window 2: AM > 10 GeV/c^ with m^o > 55 GeV/c^; 

. low neutralino mass window 3: m^o < 55 GeV/c^. 
Different selection criteria on the momentum of the tagged leptons were applied depend- 
ing on the mass window. 

Electron and muon momentum selection 

In addition to the high rejection power of the topological jet variables, lepton identification 
has been used since two opposite sign leptons of the same flavour are produced in the final 
state (see Table 2). Therefore, an electron and positron, or two muons with opposite sign 
were required, with thresholds on the momentum which depended on AM, in order to 
discriminate the selectron or the smuon pair-production signal from the SM background: 
. the momentum of the less energetic tagged lepton (electron or muon) had to be lower 

than 30 GeV/c (window 1), 70 GeV/c (windows 2 and 3); 
. the momentum of the more energetic tagged electron had to be in the intervals 
[2,40] GeV/c (window 1), [10,70] GeV/c (window 2) and [10,90] GeV/c (window 3); 
. the momentum of the more energetic tagged muon had to be in the intervals 
[2,40] GeV/c (window 1), [30,70] GeV/c (window 2) and [30,90] GeV/c (window 
3). 
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Neural network signal selection optimisation 

The following variables have been used as inputs to the neural networks: 

. the clustering variables 743, 754, ygs, computed with the CAMBRIDGE algorithm; 

. the minimum di-jet mass in the four, five and six-jet configurations; 

. the energy of the least energetic jet ■ minimum di-jet angle in four and five-jet 
configurations; 

. the thrust (only for window 3); 

. the energy of the most energetic electromagnetic cluster (only for window 3). 

The training was performed on signal samples of selectrons and smuons at centre-of- 
mass energies of 200 and 206 GeV for each analysis, and with the same statistics for two 
samples of the most important expected SM backgrounds (two and four-fermion events 
separately) . 



4.4.5 Squark search 

Searches for stop and sbottom were performed in the case of indirect decays. The eight 
quarks event topology depends strongly on AM, the mass difference between the squark 
and the Xv The same mass windows as those defined for slepton analysis (Section 4.4.4) 
were used. After the preselection and before training the neural networks, additional 
criteria were applied to select high jet multiplicity events: 

. the total multiplicity had to be lower than 40; 

. the effective centre-of-mass energy had to be greater than 0.7-y^; 

. the total electromagnetic energy had to be lower than 20 GeV (window 3 only); 

. the ln(y43) had to be greater than -6; 

. the ln(y54) had to be greater than -6.5; 

. the momentum of the less energetic tagged electron had to be lower than 16 GeV/c 
(windows 1 and 2), 20 GeV/c (window 3); 

. the momentum of the most energetic tagged electron had to be lower than 40 GeV/c; 

. the energy of the most energetic electromagnetic cluster had to be lower than 40 GeV 
(window 3 only). 
Sbottom decays produce b-quarks in the final state which may be identified with the 
impact parameter information provided by the micro-vertex detector. The event tagging 
obtained with the DELPHI algorithm for tagging events containing a b-quark [31] was 
therefore added as a sequential cut for the sbottom analysis. 

The same input variables as in the selectron and smuon searches (Section 4.4.4) were 
used in the neural network, except for the low AM analysis, where the energy of the most 
energetic electromagnetic cluster was suppressed. 

5 Results and limits 

In this section, the number of selected and expected events after the final event selec- 
tion, and the signal efficiencies obtained for each channel under study are presented. The 
results are in agreement with the SM expectation. Together with the signal efficiencies 
they were used to exclude at 95% CL possible regions of the MSSM parameter space. 
Unless otherwise stated, the limits were derived using the results from the centre-of-mass 
energies between 192 and 208 GeV. 

As already mentioned, to obtain the most conservative constraints on the MSSM 
parameter values from LLE searches, only the analyses performed considering the A133 
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coupling as the dominant one were used: in fact, if a different A^fc coupling is dominant, 
the exclusions would be at least as large as those from a dominating A133 coupling. 

After a presentation of the methods applied to derive limits, the efficiencies and the 
number of selected events are given for each channel, as well as the derived limits. 

5.1 Limit computation 

Limits on gaugino masses 

An upper limit to the number of signal events, N95, at 95% CL, was calculated according 
to the monochannel Bayesian method [32] from the number of events remaining in the 
data and those expected in the SM, summed over all centre-of-mass energies from 192 to 
208 GeV. 

The gaugino pair-production was considered for different values of tan/5 (from 0.5 to 
30), mo (between 90 GeV/c^ and 500 GeV/c^), /i (between -200 GeV/c^ and 200 GeV/c^) 
and M2 (between 5 and 400 GeV/c^); for a given set of tan/3 and mo values the (/x, M2) 
point was excluded at 95% CL if the expected number of signal, Ngxp at this point was 
greater than N95. The computations of Nexp were slightly different for LLE and IJDD 
searches, as detailed below. 

To obtain the limits on the gaugino masses with a good precision, special studies were 
performed to scan the regions of the parameter space from which the limits were deter- 
mined: the steps in M2 and /x were of 0.25 GeV/c^ and 1 GeV/c^, respectively. 

Limits on sfermion masses 

For all the sfermion searches the limits at 95% CL were derived using the modified fre- 
quentist likelihood ratio method [33]. Expected exclusion limits were obtained with the 
same algorithm where the number of observed events was set to the number of expected 
background events (absence of signal). To extract the mass limits, a branching ratio of 
100% was assumed for the i?p-conserving decay of the sfermion into a neutralino and a 
fermion. The MSSM values chosen to present the exclusion plots were tan/5 =1.5 and 
/x = -200 GeV/c2. 

The statistical errors on the efficiencies, which were between ±1% and ±3%, and on 
the expected background were used in the limit computation, for gauginos and sfermions. 
The systematic uncertainties on the signal selection efficiencies were negligeable compared 
to statistical errors in the LLE analyses. In the case of IJDD analyses, the systematic 
uncertainties on the signal efficiences were larger. Indeed, the hard gluon radiation in the 
parton shower of the Monte Carlo IJDD signal simulation is not implemented. Therefore 
this generates systematically events with background-like ymn distributions. This is the 
reason why the UDD results of the present search are conservative. 

5.2 Gaugino searches 
5.2.1 LLE scenario 

Efficiencies and selected events 

The efficiency of the selection described in Section 4.3.2 was computed from simulated 
samples at different points of the MSSM parameter space. In order to benefit from 
the high centre-of-mass energies and luminosities, all e"'"e~ -^ X^X^ and e"'"e~ -^ Xi^Xi^ 
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processes which contribute significantly have been simulated, at each MSSM point of 
this study. Then a global event selection efficiency was determined for each point. The 
efficiencies lay between 11% and 38%. 

At each selection step of the gaugino analysis, good agreement between the number of 
observed and expected background events was obtained, and no excess was observed in 
the data; at the end, 24 candidates remained in the data from 192 to 208 GeV, compared 
to 23.7 ± 0.6 expected from SM background processes (see Table 3), mainly from W+W~ 
events and the rest from other four-fermion processes. 

Limits 

The number of expected events corresponding to gaugino pair-production at each point 
of the explored MSSM parameter space was obtained by: 

i^cxp - £200 ■ Z^Ecm=192 -^Ecm ■ ^XX + ^206 " Z^Ecm=203 -^Ecm ■ ^XX 

where (t-^-^ = J2lj=i(^i^'^^~ ^ X°Xj) + o-(e+e~ -^ xtXi), ^e,^ is the integrated lu- 
minosity collected at the centre-of-mass energy Ecm, and €200 and 6206 are the global 
efficiencies determined as explained above at 200 and 206 GeV respectively. All points 
which satisfied Ncxp >N95 were excluded at 95% CL, and the corresponding excluded 
area in (/i, M2) planes obtained with the present searches are presented in Figure 7, for 
mo = 90, 500 GeV/c^ and tan/3 =1.5, 30. 

For each tan/?, the highest value of the mass of the lightest neutralino which can be 
excluded has been determined in the (/i, M2) plane for several mo values from 90 to 
500 GeV/c^; the most conservative mass limit was obtained for high mo values. The 
corresponding limit on neutralino mass as a function of tan/9 is shown in Figure 8. 
The same procedure has been applied to determine the most conservative lower limit on 
the chargino mass. The result is less dependent on tan/?, and almost reaches the kine- 
matic limit for any value of tan/5. The lower limit obtained on the neutralino mass is 
39.5 GeV/c^, and the one on the chargino mass is 103.0 GeV/c^. 



5.2.2 UDD scenario 

Efficiencies and selected events 

At the end of the analysis procedure to search for gauginos described in Sections 4.4.2 
and 4.4.3, no significant excess of events was seen in the data with respect to the SM 
expectations. Figure 9 shows the number of expected events from the SM and the number 
of observed events as a function of the average signal efficiency obtained with all simulated 
masses for the N3 and C2 analyses after a step- by-step cut on the neural network output. 
For neutralino pair-production, the efficiencies were typically around 30-60% at the 
values of the optimized neural network outputs, depending on the simulated masses. For 
chargino pair-production, the signal efficiencies were between 10% and 70%. The ex- 
pected and observed numbers of events for both analyses are reported in Table 3 for each 
mass window. 

Limits 

The signal efficiency for any values of the Xi and x^ masses was interpolated using an 
efficiency grid determined with signal samples produced with the full DELPHI detec- 
tor simulation. The number of expected events Nexp has been computed separatly for 
neutralino and chargino pair productions. 
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1\T ^ Y~^I^cm— 202 r> ^ \ r- Y~^I^cm — 208 r> 

INexp - £200 ■ Z^Ecm=192 -^Ecm " ^^XX + ^206 " Z^Ecm=203 ^^Ecm " (^XX 

where a^-^ = o"(e"'"e~ -^ XiXi) or ^(e+e^ -^ Xi^xr)) ^.nd £200 and 6206 are taken from the 
efficiency grids. 

Using the six-jet and the ten-jet analysis resuhs, an exclusion contours in the (/i, M2) 
plane at 95% CL were derived for different values of mo (90 and 300 GeV/c^) and tan/? 
(1.5 and 30), as shown in Figure 10. In the exclusion plots the main contribution comes 
from the study of the chargino indirect decays with the ten-jet analysis, due to the high 
cross-section. The six-jet analysis becomes crucial in the exclusion plot for low tan/?, low 
mo and negative /i. The lower limits on the mass of the lightest neutralino and chargino 
are obtained from the scan on tan/5 from 0.5 to 30. The lower limit on the neutralino 
mass of 38.0 GeV/c^ is obtained for tan/? = 1 and mo = 500 GeV/c^ (Figure 11). The 
chargino is mainly excluded up to the kinematic limit at 102.5 GeV/c^. 

5.3 Slepton searches 
5.3.1 LLE scenario 

Efficiencies and selected events 

The efficiencies of the slepton analysis described in Section 4.3.3 were between 18% and 
38% for the sneutrino indirect decay channel, depending only on the neutralino mass. 
The efficiencies were higher for the final states obtained in indirect decay of charged 
slepton pairs, due to the presence of two additional leptons. They ranged from ~ 20% 
(m^^o = 15 GeV/c^) to 43% for stau pairs; they were of the same order but up to ~ 5% 
higher for selectron pairs, and ranged from ~ 25% (m^^o = 15 GeV/c^) to 64% for smuon 
pairs. For the direct decay of u^, the analysis efficiencies lay in the range 27-36%, de- 
pending on the sneutrino mass. 

At the end of the selection, 11 events remained in the data compared to 8.1 ± 0.3 
expected from the SM processes (Table 3). The background was mainly composed of 
four-fermion events, in particular from W pair-production. 

The efficiencies of the z/,-//,- direct decay analysis (see Section 4.3.4) varied with the 
i)r mass and ranged from 45% to 51%. At the end of the selection, 6 candidates were 
obtained for 6.3 ± 0.4 expected (see Table 3). 

Limits 

To derive limits on slepton masses, the results of the search described above were com- 
bined with those obtained with data at y/s = 189 GeV [3]. 

For charged slepton indirect decay, the areas excluded in the m^o versus m^ planes 
are plotted in Figure 12. 

As was explained in Section 2.3, a pair of selectrons can be produced in the t-channel 
via neutralino exchange. With the MSSM parameters fixed to derive limits, the e"'"e~ 
cross-section is higher than the fx^jl~ and f^f~ ones. So, though the analysis effi- 
ciencies for the smuon pair-production were higher, the excluded area in case of the 
selectron pair search is the largest; the smallest is obtained for the f+f^ production. For 
AM > 5 GeV/c^ the limits on the slepton mass are 94 GeV/c^, 87 GeV/c^ and 86 GeV/c^ 
for the e, fx and f, respectively, and become 95 GeV/c^, 90 GeV/c^ and 90 GeV/c^ if the 
neutralino mass limit is taken into account. 

The results of the search for the indirect decay of the z/ were used to exclude areas 
in the vHyXt versus m^ planes, as shown in Figure 13. These exclusion areas are also 
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valid for all the Ajj^ couplings. As already mentioned in Section 2.3, the z/c^c cross- 
section can be enhanced compared to the z/^z/^ and z/t-z/,- cross-sections if production via 
a chargino exchange is possible and the excluded area depends on the chargino mass. For 
AM > 5 GeV/c^, the limit on the z/^ and z/^ mass is 82 GeV /(?, and is 85 GeV/c^ if the 
limit on the neutralino mass is taken into account. These limits are 96 and 98 GeV/c^ 
respectively for u^. 

The results of the searches for 4r and 2e2r final states, from sneutrino pair direct 
decays, were combined to obtain lower limits on the sneutrino mass. The results from 
the 4r search were used to derive limits on z/e and on z/^, those from the 2e2r to derive 
limits on z/,-. The limits obtained are respectively 96 GeV/c^, 83 GeV/c^ and 91 GeV/c^. 

5.3.2 tJDD scenario 

Efficiencies and selected events 

No significant excess has been observed in the output signal node distributions for any 
analyses. The signal output distributions for the selectrons and smuons analyses are 
shown in Figure 14 for the medium AM analyses. The signal efficiency of the slepton 
analyses described in Section 4.4.4 was evaluated at each of the simulated points for the 
two centre-of-mass energies (200 and 206 GeV). Efficiencies for the signal (selectron and 
smuon) were in the range from 5-40%, for small mass differences and small neutralino 
mass, and increased up to 60% for medium AM analyses. The 5% efficiency was obtained 
for the AM = 5 GeV and for a neutralino mass of 45 GeV. This efficiency increased rapidly 
with AM and with the neutralino mass. 

No excess of data with respect to the SM expectations was observed for the selec- 
tron and smuon analyses; the numbers of events observed and expected from background 
contributions are shown in Table 3. The remaining background comes mainly from four- 
fermion processes. 

Limits 

From the selectron and smuon pair-production searches exclusion domains have been 
computed in the m^o versus m^ plane (Figures 15). The M2 value was fixed for each 
neutralino mass. For AM> 5 GeV/c^, the lower limit on the right-handed selectron mass 
was 92 GeV/c^, and the lower limit obtained for the right-handed smuon was 85 GeV/c^. 

5.4 Squark searches 
5.4.1 LLE scenario 

Efficiencies and selected events 

The selection efficiencies of the analysis described in 4.3.5 varied with the stop mass and 
with the mass difference between the stop and the lightest neutralino. They lay around 
30% for most of cases, and around 18-20% for low neutralino masses. 

After the selection procedure, 35 events remained with 35.4 ± 0.6 expected from back- 
ground contributions, mostly coming from W"'"W~ production. 

Limits 

From the study of the stop indirect decay, a lower limit on the stop pair-production 
cross-section was derived as a function of the stop and neutralino masses. Considering 
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two cases of stop mixing (no mixing and mixing angle = 56°), the exclusion limit was 
derived in the m^o versus m^ plane, as shown in Figure 16. With no mixing, the lower 
bound on the stop mass is 88 GeV/c^, valid for AM > 5 GeV/c^. If the mixing angle 
is 56°, the lower bound on the stop mass is 81 GeV/c^, for AM > 5 GeV/c^, and becomes 
87 GeV/c^, taking into account the lower limit on the mass of the lightest neutralino. 

5.4.2 tJDD scenario 

Efficiencies and selected events 

The final selection of candidate events was based on the signal output values of the 
neural networks for the stop and for the sbottom. The cut on the neural network variable 
has been relaxed for sbottom due to the effect of the b-tagging selection. The signal 
outputs of the neural network for the multi-jet stop and sbottom analyses (window 3) 
are shown in Figure 17. The signal efficiencies of the neural network analyses described 
in 4.4.5 were evaluated at each of the evenly distributed simulated points in the plane of 
stop (sbottom) and neutralino masses. Efficiencies for the signal after the final selection 
were in the range from 10-20%, for small or large mass differences between squark and 
neutralino, up to around 50% for medium mass differences. 

The numbers of events observed and expected from SM processes are shown in Table 3. 

Limits 

The resulting exclusion contours for stop and sbottom can be seen in Figures 18 and 19. 
The lower limit on the mass of the left-handed stop, assuming the neutralino mass limit 
of 38 GeV/c^, is 87 GeV/c^ for AM > 5 GeV/c^. The lower limit on the stop mass with 
the mixing angle = 56° was 77 GeV/c^ under the same assumptions. The lower limit 
obtained for the left-handed sbottom assuming a neutralino mass limit of 38 GeV/c^, 
was 78 GeV/c^ for AM > 5 GeV/c^. The sbottom pair-production cross-section in the 
case of the mixing angle corresponding to maximum decoupling from the Z boson (68° 
for sbottom) was too low to cover a significant region of the mass plane by the excluded 
cross-section of the sbottom analyses. 

6 Summary 

A large number of different searches for supersymmetric particles with the assumption 
of Rp violation via AjjfcLjLjEfc or Af -^UjDjDfc terms have been performed on the data 
recorded in 1999 and 2000 by the DELPHI experiment, at centre-of-mass energies between 
192 and 208 GeV. No significant excess has been observed in any of the channels. Limits 
on the pair-production of sparticles have been derived. These limits were converted 
into limits on sparticle masses and excluded regions in the MSSM parameter space. Mass 
limits are summarized in Table 4, together with the assumptions under which these limits 
are valid. 
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v^ (GeV) 


192 


196 


200 


202 


<204.9> 


<206.6> 


<206.6> 


C (pb-i) 


25.1 


76.0 


83.3 


42.5 


73.7 


85.4 


51.8 



Table 1: Data collected by DELPHI in 1999 and 2000: the integrated luminosities corre- 
spond to the data actually used in the present analyses after the run selection. The last 
column refers to the integrated luminosity collected with one sector of the TPC off. 
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LLE: multi-lepton topologies 



analysis name 


final states 


direct 
decays of 


indirect 
decays of 


gaugino 


2r + n£ + mj + ^ 
(n>2) (m>0) 


9090 

AlAl 


xtxi 
xU] 


slepton 


2r + 2£ + ^ + p£ 
(p = or 2) 

TTTT 


l^ei^e 


z>^, £+£- 


sneutrino tau 


eeTT 


VrVr 




squark 


2t + 21 + $+ 2j 




ql 


UDD: multi-jet topologies 


analysis name 


final states 


direct 
decays of 


indirect 
decays of 


neutralino 


6j 


x?x?, xtx. 




chargino 


lOj 




xXxl 


slepton 


6j + 2i 




I+I- 


squark 


8j 




ql 



Table 2: The multi-lepton and multi-jet visible final states which correspond to the analy- 
ses described in this paper, when one LLE or UDD term is dominant. The corresponding 
pairs of produced sparticles that may have given rise to them are indicated. For the 
LLE cases, only topologies produced with decays via A133 are considered (see text), and 
for the UDD cases, i = e or fi. 
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LLE: multi-lepton topologies 



Analysis 






192- 


-202 GeV 


203- 


-208 GeV 




name 






observed 


expected 


observed 


expected 


gaugino 






15 


13.4 


± 0.4 


9 


10.3 


± 0.4 


slepton 






4 


4.2 


± 0.2 


7 


3.9 


± 0.2 


sneutrino r 






3 


3.5 


±0.3 


3 


2.8 


± 0.2 


squark 






19 


19.9 


± 0.5 


16 


15.5 


± 0.4 






UDD: multi-jet topolo, 


gies 








Analysis 


Mass windows 




192- 


-202 GeV 


203- 


-208 GeV 




name 


(GeV/c2) 




observed 


expected 


observed 


expected 


neutralino 


15 < m^ < 45 




134 


126.0 


± 13.0 


121 


119.3 


± 8.8 




45 < m^ < 75 




192 


172.5 


± 8.2 


167 


164.7 


± 5.7 




75 < m^ 




97 


103.3 


±3.6 


82 


91.7 


± 2.3 


chargino 


AM < 10 




187 


181.1 


± 5.9 


156 


171.7 


± 5.2 




AM > 10 




22 


25.6 


± 1.1 


20 


23.5 


± 1.0 


slepton 


m^ > 55, AM < 10 


9 


5.6 


± 0.2 


1 


6.2 


± 0.2 


(g) 


m^ > 55, AM > 


10 


1 


2.0 


± 0.1 


5 


2.3 


± 0.1 




15 < m^ < 55 




1 


1.6 


± 0.1 





1.8 


± 0.1 


slepton 


m^ > 55, AM < 10 


7 


5.7 


± 0.2 


5 


6.4 


± 0.2 


il^) 


m^ > 55, AM > 


10 


4 


3.3 


± 0.2 


1 


3.5 


± 0.2 




15 < m^ < 55 




2 


2.0 


± 0.1 


1 


2.3 


± 0.1 


squark 


m^ > 55, AM < 10 


42 


39.4 


± 0.6 


38 


40.4 


± 0.6 


(t) 


m^ > 55, AM > 


10 


13 


10.1 


±0.3 


8 


9.5 


±0.3 




15 < m^ < 55 




30 


26.3 


±0.5 


25 


25.2 


±0.5 


squark 


m^ > 55, AM < 10 


10 


11.9 


± 0.4 


13 


12.0 


± 0.4 


(b) 


m^ > 55, AM > 


10 


4 


5.4 


± 0.2 


7 


4.8 


± 0.2 




15 < m^ < 55 




6 


4.6 


± 0.2 


8 


4.2 


± 0.2 



Table 3: Observed and expected numbers of events for all LLE and UDD analyses. 
Although not explicitely written, AM is always greater than 5 GeV/c^. 
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SUSY 
particle 


Comnients about validity conditions 


Mass limii 
LLE 


t (GeV/c^) 
UDD 


X? 

xt 


Validity conditions for gauginos: 

90 < mo <500 GeV/c2, 0.7 < tan/3 <30 , 

-200 < ;u < 200 GeV/c^ and <M2 < 400 GeV/c^ 


39.5 
103.0 


38.0 
102.5 




Validity conditions for sfermions: 
/i = -200 GeV/c2 and tan/3 = 1.5, 
BR(f^ f X? )=1, AM>5 GeV/c2 






gr 


x5 niass limit not used 


94 
95 


92 




X? mass limit used 


92 


m 


X? mass limit not used 


87 
90 


85 




X*]* mass limit used 


87 


tr 


Xi mass limit not used 


86 
90 


- 




x5 mass limit used 


- 




X? mass limit not used 


96 
98 
96 


— 


^e 


X? mass limit used 


- 




direct decay only 


X 




Xi mass limit not used 


82 
85 
83 


- 


^^ 


X? mass limit used 


- 




direct decay only 


X 




X? mass limit not used 


82 
85 
91 


- 


Vr 


X? mass limit used 


— 




direct decay only 


X 




Xi mass limit not used, no mixing 


88 
92 

81 

87 


81 


t 


Xi mass limit used, no mixing 


87 




X? mass limit not used, minimal mixing 


67 




X? mass limit used, minimal mixing 


77 


b 


X? mass limit not used, no mixing 
X? mass limit used, no mixing 





78 
78 



Table 4: Sparticle mass limits at 95% CL from the DELPHI R^ violation pair- 
searches of supersymmetric particles; x : the decay channel is not possible; - 
channel is not covered. These results are valid if m^o > 15 GeV/c^, except for 
limits derived from the v direct decay which does not depend on m^o . 



production 

the decay 

the V mass 



27 




r± {t+) 
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b) 
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Figure 1: Diagrams of sparticle direct decays, a: slepton direct decay via LLE term; 
b: squark direct decay via IJDD term; c: neutralino/chargino direct decay via any 
i?p-violating trilinear term. 








f f) 




Figure 2: Diagrams of sparticle i?p-conserving decays. To get the whole chain of the 
sparticle indirect decay, the LSP (lightest neutralino) has to undergo a direct i?p-violating 
decay, a, b: slepton decays into a lepton and the lightest neutralino; c: squark decay into 
a quark and the lightest neutralino, the hatched disk means decay beyond tree level for 
the stop case; d, e, f: examples of gaugino decays. 
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Figure 3: LLE: sparticle pair search in multi-lepton channels - Event variable (see text) 
distributions before the criteria applied on the acollinearity (preselection) (a) and af- 
ter (b-f) the preselection. The simulated signal corresponds to x'Jx'j' production, with 
m^o = 90 GeV/c^ (the normalisation is arbitrary). 
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Figure 4: LLE: sparticle pair search in multi-lepton channels - Distributions of variables 
used in the four analyses, as described in the text: a) gaugino, b) slepton, c) sneutrino 
tau, d) squark. The simulated signals correspond to production of a) X2X1 and xtXi 
(m^o = 50 GeV/c^, m^o = 85 GeV/c^, m~± = 95 GeV/c^), b) u^3^ (m^^ = 100 GeV/c^), c) 

95 GeV/c^), d) tt (m^ = 95 GeV/c^) (the signal normalisation is arbitrary). 
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Figure 5: UDD: the total energy associated to charged particles (left) and -ln(y43) (right) 
distributions in the multi-jet slepton analyses at the hadronic preselection level. 
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Figure 6: UDD: the -ln(y54) (left) and -ln(y65) (right) distributions in the multi-jet 
sfermion analyses at preselection level. 
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Figure 7: LLE: regions in fi, M2 parameter space excluded at 95% CL by the neutralino 
and chargino searches for two values of tan/? and two values of mo. 
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Figure 8: LLE: excluded lightest neutralino mass as a function of tan/5 at 95% CL. This 
limit is independent of the choice of the generation indices i,j,k of the Xijk coupling and 
is for values of mo between 90 and 500 GeV/c^. 



202 - 208 GeV 



202 - 208 GeV 




80 100 

Efficiency (%) 



10- 



10' 



10 




10 20 30 40 50 60 70 80 

Efficiency (%) 



Figure 9: UDD: number of expected events (continuous line) and data events (black 
dots) versus average signal efficiency for the high neutralino mass search N3 (plot a) and 
for the large AM chargino search C2 (plot b) for centre-of-mass energies between 202 and 
208 GeV. 
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Figure 10: UDD: regions in /i, M2 parameter space excluded at 95% CL by the neutralino 
and chargino searches for two values of tan/3 and two values of mo. 



34 



(N 



70 p- 

65 r- 



DELPHI 4 = 




Figure 11: UDD: lightest neutralino mass excluded at 95% CL as a function of tan/3. 
This limit was obtained for mo= 500 GeV/c^. 
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Figure 12: LLE: excluded regions at 95% CL in the m-j^o versus m^^ planes for iji 

pair-production, with BR(£fl ~^ X? ^) = 100%, and neutralino decay into leptons. The 
plots show the exclusion (filled area) for selectron, smuon and stau. The black contour 
is the corresponding expected exclusion at 95% CL. 
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Figure 13: LLE: excluded regions at 95% CL in the ni-^o versus vny planes for v^ (left) 
and Vp^^T (right) pair-production, with BR(i> — ^ Xi '^) = 100%, and neutralino decay into 
leptons. The black contour is the corresponding expected exclusion at 95% CL. 
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Figure 14: UDD: the neural network signal output distributions for the selectron (left) 
and smuon (right) window 2 analyses. The cuts on the neural network output variable 
were chosen for the final selection at 0.83 (selectrons) and 0.92 (smuons). 
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Figure 15: UDD: excluded domains at 95% CL in the m^o versus m^ planes for selectron 

(left) and for smuon (right) pair-production with BR(£ij —>■ Xi ^) = 100% and neutralino 
decay into jets (filled area). The superimposed contours show the expected exclusion at 
95% CL. 
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Figure 16: LLE: exclusion domains at 95% CL in the m^o versus m^ plane for the stop 
pair-production, with BR(t — *> cx°) = 100% and neutralino decay into leptons. The 
plots show the exclusion (filled area) for the lightest stop for no mixing (left) and for the 
mixing leading to the maximal decoupling to the Z boson (right). The black contour is 
the corresponding expected exclusion at 95% CL. 
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Figure 17: UDD: the neural network signal output distributions for the multi-jet stop 
(left) and sbottom (right) window 3 analyses. The cuts on the neural network output 
variable were chosen for the final selection at 0.86 and 0.75 for stop and sbottom respec- 
tively. 
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Figure 18: UDD: exclusion domains at 95% CL in the m^o versus m^ plane for the stop 
pair-production, with BR(t — * ex?) = 100% and neutralino decay into jets. The plots 
show the exclusion (filled area) for the lightest stop for no mixing (left) and for the 
mixing leading to the maximal decoupling to the Z boson (right). The black contour is 
the corresponding expected exclusion at 95% CL. 
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Figure 19: UDD: exclusion contours at 95% CL in the m^o versus nij;, plane for the 

sbottom pair-production, with BR(b -^ ^Xi) = 100% and neutralino decay into jets. 
The plot shows the obtained exclusion (filled area) for the lightest sbottom in the case of 
no mixing. The black contour is the expected exclusion at 95% CL. 



